EIGENVALUE PROBLEM IN TWO DIMENSIONS FOR AN IRREGULAR 

BOUNDARY 
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An analytical perturbative method is suggested for determining the eigenvalues of the Helmholtz 
equation (V 2 + k 2 )tp — in two dimensions where ip vanishes on an irregular closed curve. We can 
thus find the energy levels of a quantum mechanical particle confined in an infinitely deep potential 
well in two dimensions having an irregular boundary or the vibration frequencies of a membrane 
whose edge is an irregular closed curve. The method is tested by calculating the energy levels for 
an elliptical and a supercircular boundary and comparing with the results obtained numerically. 
Further, the phenomenon of level crossing due to shape variation is also discussed. 



I. INTRODUCTION 

The energy levels of a quantum particle confined in a 
2D regular box can be solved exactly only in the cases of 
a square and a triangle and in the limiting case of a cir- 
cle. While the determination of the energy levels for the 
circular or the square boundary is a trivial exercise, the 
problem of the triangular boundary is more formidable 
The corresponding problems in the classical regime 
can be the flow of liquid through a pipe of polygonal 
cross-section or the free vibration of a membrane (with a 
fixed boundary) of polygonal shape. The classical prob- 
lems, like their quantum counterparts, are amenable to 
simple analytical treatments only in the cases of a cir- 
cle, a square and a triangle. The problem of a regular 
polygonal box has been solved by perturbing about the 
equivalent circle and the results have been quite accu- 
rate @. The same problem has been solved by Cureton 
and Kuttler [3[ in the context of vibration of membranes. 
Here we address the problem of finding out the energy 
eigenvalues when the boundary has no simple geomet- 
ric shape. The Schroedinger equation for a particle of 
mass m and energy E confined in an infinitely deep 2D 
potential well is 



h 2 2 

— V 2 ^ = E$ 
Ira 



which can be recast as 



(V 2 + k 2 )i> = 



(la) 



(lb) 
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where k — y ^jpr ■ Thus the problem boils down to 
solving the Helmholtz equation with the Dirichlet condi- 
tion ip = on the 'irregular' boundary. Exact solutions 
can be obtained only in a few special cases. The standard 
procedure is to choose a curvilinear coordinate system 
suited to the geometry of the problem and employ the 
method of separation of variables. For a boundary having 
an irregular shape no particular coordinate system will be 
useful. Hence, we resort to perturbative methods to solve 
the problem. Here we will perturb the boundary about 
a circle so that in our problem solutions can be obtained 
in the form of corrections to the solutions for the circular 
boundary. Till now, most of the effforts at finding out 
the eigenvalues of the Helmholtz equation for an irregu- 
lar boundary have been numerical. Mazumdar 0, 0, E| 
reviews the approximate methods invoked for this prob- 
lem. In addition to the extensive summary of theoretical 
results, Kuttler and Sigillito [?J also gives a comprehen- 
sive review of the different numerical methods employed. 
More recently, Amore 0] gives a numerical recipe using 
a collocation approach based on little sine functions. As 
far as analytical works are concerned, Rayleigh [9( and 
also Fetter and Walecka jTOj] finds the ground state energy 
eigenvalues for a vibrating membrane. A general formal- 
ism has been suggested by Morse and Feshbach [ll| using 
the Green functions. Parker and Mote [l2j have put for- 
ward a perturbative method for finding the eigenvalues 
and the eigenfunctions through fifth order. A similar 
method has been proposed by Nayfeh[l3]. However, the 
eigenvalues are found out only to the first order. Read 
[1J| has also suggested ageneral analytical approach to 
the problem. Bera et al [15l | have proposed a perturbative 
approach to the problem but failed to express the solu- 
tions in a closed form. Our approach is similar in spirit 
to that of Bera. Here we present a closed form solution 
to the problem in a more systematic and efficient man- 
ner. The method is tested by comparing the analytical 
results with those obtained numerically for a supercircu- 
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lar and an elliptical boundary. Further, the phenomenon 
of energy level crossing as induced by the shape variation 
is also dealt with for both the boundaries. In section II 
we set up our general scheme and in section III we apply 
it to the case of a supercircle and an ellipse. A short 
conclusion is presented in section IV. 



II. PERTURBATION ABOUT THE 
EQUIVALENT CIRCLE 

Given any r(9) defining the boundary in 2D enclosing 
an area A we first construct a circle of radius Ro such 
that 



Now, as a first approximation, the energy E a of 
the particle confined by r{9) will be that of a particle 
enclosed in a circle of radius R a 



E n = 



2„2 
tn.n 



h 2 p. 



2mRl 



(6) 



with p mn = k m n R being the nth node of the mth 
order Bessel function. The next step is to improve 
upon the 'equal area' approximation by perturbing the 
equivalent circle and finding out the first and the second 
order corrections to the eigenvalues. 

We expand ip and E as 



A = irRl 



(2) 



We can then expand r(9) about Rq in terms of Fourier 
series at different orders of smallness (denoted by A) as, 



r(9) = Ro 



(3) 



where, 



/( <J )(6») = YVC^' cosn<9 + S { n a) sinn<9) 



(4) 



Here, for simplicity, we have considered only a one pa- 
rameter dependence of r(9), which thus represents a fam- 
ily of curves which reduce to the equation for a circle in 
the limiting case A — > 0. A should be much smaller than 
unity which ensures that the variation of r{6) with 9 is 
small enough to permit the use of perturbativc methods. 
If we now calculate the area using (3), (4) and equate it 
with (2) then the following constraint relations among 
the Fourier coefficients are arrived at, 



OO u — 1 

E E [ck v) c£- v) + 2<t ) cf + spsw 



ip = ipo + Xipi + A ip2 + 



E = E + XE 1 + X 2 E 2 + 



(7a) 



(7b) 



Using (|7aj) . ((7b]) in (Tlaj) , equating the coefficients of 
different powers of A to and after some rearrangement 
we arrive at the set of equations, 



2mE n 



h 2 



-ipo 



(V 2 



(V 2 + kl)^ = - 



2m 



2mEi 



(8a) 



(8b) 



■ = -jp{Eii>i + £ 2 ^ ) (8c) 



jaj) can readily be identified as the equation for the cir- 
cular boundary with as the eigenfunction correspond- 
ing to energy Eq. 
The boundary condition is 

V(i?o + Ai?o/ (1) + A 2 i? / (2) + -..) = 



Taylor expansion about r = Rq, use of (|7ap and equating 
the coefficients of different powers of A to yields 



AC, 



(5a) 



^0(^0) = 



(9a) 



In particular we have, 



MRo) + R f {1) %(Ro) = 



(9b) 



C« = 



(5b) 



4Gf = -El^ 132 + ^ 



> 2 i 



(5c) 



+ i? / (2) V^(i?o) + ^o/ (1) Vo (Ro) = 
We discuss separately the cases I = and 1^0. 



(9c) 
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A. Calculation of Energy for / = State 

For the 1=0 state, 

i>o = NJo(p) 



(10) 



where Jo is the 0th order Bessel function, and N is 
the normalisation constant. E is obtained from ^ with 
m = 0, and an appropriate n. 

The first order correction to the wave function, obtained 
as a solution to (I8bl) is, 



ipi = /(a p cos p9 + a p sin p9)J p + a J - ttttNJi (11) 

, 2E 

p=i 

where the last term is the particular integral to (I8bl) . 
Incorporating pip in (|9b[) and separately matching the 
coefficients of the cosine and the sine terms we have 



-Po, 



Jp(P0,n) 



-Po 



P Jp(P0,n) 



(12a) 



(12b) 



/ ^^p+fc^p + °|fe-p| L/ P ^k-p^p + °p-fc D p \ 

(14b) 



= -po, n Jo[NS^+a S ( k 1) }+ 



PQ.uJqN ( 1 P0,n J'p 



p=l 

(14c) 

As before, the remaining constant bo can be deter- 
mined by normalising the wavefunction upto the order 
of A 2 . 



B. Calculation of Energy for I 7^ state 

The I ^ states come in 2 varieties, 



"•• - v '-'"">l Zw 



For simplicity, we assume that S„ = for all a. We 
will work with 



(15) 



E 1 = 



(12c) 



The remaining constant ao can be found out by nor- 
malising the corrected wavefunction over the enclosed 
area. However, that is not required right now for our 
purpose. (12c) implies that there cannot be any correc- 
tion to the energy in the first order. So any possible 
correction to the energy can only come from the second 
or higher orders. 

In a similar fashion the second correction to the wave 
function as a solution to (l8c|) with E\ = is found out 
to be, 



pE 2 

2E V; 



which, when introduced in (|9c[l . now yields 



p cos p9 + bp smp9)J p + bo Jo 

P =i 



iVJi (13) 



Eo — En 



Lp=i 



PO.n 



Jp(P0,n) 1 
Jp(P0,n) 2 



(2) 



2C, 



(14a) 



ipo = NiJi(p) cos 19 



(16) 



The result for the other case will be similar. The first 
correction to the wavefunction obtained as a solution to 

ia is, 



ipi = cipJp cos p9 + I aiJi 

p=0,p#i 



Eip 

En 2 



NiJi+i cos 19 



Following a similar procedure as that for the ground 
state we now have, 



Ex = -C$>Eo 



Pi ,3 AT j i(pi,j) tnW j-M 1 ) 



(18a) 



- N IJ^)K+i+ C \p-1\) forp^0,Z (18b) 



„ _ Phi Ar J l(Pu) r (i) 
2 MPij) 



(18c) 



with ai being obtained from the normalisation condi- 
tion. The second order corrections yield, 



b k = -po,nJo[NC i *>+a Ci 1 >}+ 



P0,nJ p 

2J<n 



Lt,o 



cos m9 



4 



r (i)2P T E 2 



En 



Ni^ cos W 



(19a) 



III. APPLICATION TO SIMPLE CASES 
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^ n — 1 



-4 2) + E 

7i— l,n^l 



(1) j^M \2 Pl,3 J 'n{Pl,i) i c (l)2 Pi, J Jq(PI,j ) 



^ ' " + ' ' " ''' 2J n ( P/J ) ^ J (p ;j ) 

(19b) 

The constants 6 m can also be determined as in the case 
of the ground state. However, that is not needed for our 
present purpose. For the case, 



i) Q = NiJi{p) sinZ6> 
similar calculations result in, 



(20) 



The general formalism having been outlined above we 
now estimate the energy levels of a supercircle and an el- 
lipse where direct comparison with the numerical results 
can be made. 

B. Particle Enclosed in a Super circular Enclosure 

Super circle is a special case of a super ellipse [l6[ 
whose equation is given by the Lame equation, 



a" 



1 



(22) 



with n > and rational, a and b being positive real 
numbers. They are also known as Lame curves or Lame 
ovals [l7| • Super ellipses can be parametrically described 
as, 



ipi= E CL p JpS\np9 + ( aiJi - ^-^N t Ji +1 ) sinW 



--i, V jkl 



E x = C^Eq 



(21a) 



(21b) 



^ = ^ Nl f^^- C lP-^ fOT ^ M (21C) 



x = a cos™ (t) 



y = b sin" (t) 



(23) 



(24) 



Different values of n gives us closed curves of different 
shapes. For n > 2 we consider only the real positive val- 
ues of cos " (t) and sin™ (t) for ^ t ^ | and use the sym- 
metry of the figure to continue to the other quadrants. 
We are interested in the case a = b which corresponds to 
a supercircle. In polar coordinates the equation for the 
supercircle is, 



a = 



(21d) 



(cos™ 9 + sin™ 0) » 
and the radius of the equal area circle is, 



(25) 



sin rat 



+ 



Ni-svaie 



(21e) 



7-i ^(1)2 , oo 

^ 2 _ L '2i | 1 y r (i) ( g r (D r (i) _ r U) \_2C {2) 



c. 



(2) 



(1) ^(1) A 2 Pl,jJn(Pl,j) 
+1 U \n-l 



n—O.n^l 



(21f) 



Rn = a 



2 [r(i) 



[r(; 



(26) 



With the deformation parameter <5 = 2 — n one arrives 
at after some rigorous calculation, 



r = Ro 



1 + C in cos 4n6» + C in cos 4n6> 

71 — 1 71—0 

where the Fourier coefficients are found to be 



(27) 



1 



4n(4n 2 - 1) 



32 \_ 8 9 / ' 



cf } is given by (5c) to be .0070205 and cf ) 



0357983. 
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Using these Fourier coefficients, the first few energy lev- 
els are calculated for the supercircular boundary in the 
range — 1 ^ 8 ^ 1, and compared with the numerically 
obtained values. This is shown in FIGQ] 

where the numerical results are shown by discrete 
points and the analytical ones by the continuous lines. 
The fact that even for such a wide range of 8 the 
analytical results are in fairly good agreement with 
those obtained numerically does indeed justify the 
validity of our formalism. Furthermore, it is to be noted 
that the energy level corresponding to the unperturbed 
wavefunction i/jq = N2J2COS29 is strongly affected com- 
pared to the others and crosses over to its counterpart 
tpo = N2J2 sin 29 at 8 — 0. This crossing of energy levels 
is solely induced by the variation in the shape of the 
boundary of the potential well. 

C. Particle Enclosed in an Elliptical Enclosure 

The determination of the eigenvalues of the Helmholtz 
operator in 2D with an elliptical boundary has been in- 
vestigated extensively. However, most of the efforts have 
been directed at the numerical estimation of the eigen- 
values d, H3| ■ An analytical method has been sug- 
gested by Wu and Shivakumar [2lJ. Here we propose a 
simpler approach to the problem. The equation for an 
ellipse with semiaxes a and b, in polar coordinates is, 



(28) 



'1 - (1 - £) cos 2 e 

Defining the deformation parameter , 



A = 



a — b 
a + b 



(29) 



28]) can be recast as, 



r = R [1 + X cos 29 - i A 2 + ^ A 2 cos 40 + A 3 () + . ..] (30) 



Comparing with our general Fourier series we ob- 
serve that, dp = 1, dp = -\, dp = §. 

(fT^^ . (fT^a^ . (fTSalt . (fT9Bh . (f^TBb . (121^ on using (5c) then 
gives, 



£1 



Ji cos 19 
Ji sin 19 



E 2 



Ji cos 19 
J 1 sin 19 



_ 1 1 Pl,jJg(Pl,j) 
2 Jb(pu) 



8_a 
2 



E0S21 



E Q 8i 



2J P (pi,j) 



\P-1\ 



(31a) 



(31b) 



The results for the elliptical boundary are shown in 
FIGUl From FIGO it is seen that as in the case of 
the supercircle here also the J\ cos 9 state is strongly 
affected by the boundary perturbation and crosses 
over to its counterpart Jisin# at 8 = 0. However, 
quite interestingly, the J2 states do not cross but are 
rather repelled by each other. They touch each other 
tangentially at 8 = 0. While for one of these states the 
analytical method works quite well it has a restricted 
validity for the other one. In fact, we compared the 
energy levels for the first 10 states and found out the 
agreement between the analytical and the numerical 
results to be quite satisfactory except when the levels 
repels each other. This phenomenon of level repulsion, 
also goes by the name of "loci veering" . In case of 
level repulsion the validity of the perturbation theory is 
restricted to a small range in 8 (e.g. \8\ ^ 0.08). This is 
in sharp contrast to the case where there is no repulsion 
in which case the agreement with the perturbation 
theory persists over a wide range. 

IV. CONCLUSION 

One of the principle virtues of the method proposed 
is its generality. With slight modification, the formalism 
can readily be adopted to study the shape dependence of 
the eigenvalues of a vibrating membrane with Dirichlet 
conditions on an irregular boundary. The approach can 
also be useful in studying the modes of propagation of 
electromagnetic waves in waveguides with irregular cross 
sections. Another potential area where this formalism 
might be useful is in the study of quantum dots. This 
field has been an area of vigorous research for the past 
few years. 2D quantum dots are generally taken to have 
a circular symmetry. However, in practice such a sym- 
metry can not be strictly ensured. There is bound to be 
small deviations from exact circular symmetry. Hence, 
probes have been constructed to investigate the shape 
of the dots [22|, HH, [24|. As shown in this paper the 
energy eigenvalues of a particle confined in 2D in an in- 
finitely deep potential well will essentially depend upon 
the shape of the confining region. Hence a study of the 
shape dependence of the energy levels might proove to 
be useful in shedding light upon the actual shapes of the 
dots. Another significant aspect of our formalism is the 
use of the general Fourier series to express the deviation 
of the boundary from a circular one which allows us to 
treat any sort of boundary within the limit of small per- 
turbation for which our formalism is valid. Further, the 
solutions are found out exactly in a closed form without 
any major approximations which is indeed remarkable. 
The case of the supercircular membrane shows that even 
for quite large perturbations the method yields satisfac- 
tory results. The accuracy of the method can be still im- 
proved by including higher order corrections. In fact, we 
have also found out the third order corrections, although 
the results are not included here. On the contrary, the 
case of the elliptical boundary points out to the failure of 
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the perturbation theory whenever the energy levels exibit 
loci veering. This provides potential topics for future in- 
vestigations. Another point which we want to emphasis 
here is that the success (and also the efficiency) of the 
formalism depends to a large extent upon the judicious 
choice of the deformation parameter A. For the case of 
the ellipse we defined A to be A — whereas the ec- 
centricity e would seem a more appropriate candidate for 
A. For the elliptical boundary we have considered defor- 
mations upto the extent where a : b = 2 : 1 for which 
A = 0.333. Had we formulated the problem in terms of 
the eccentricity the same deformation would have led to 
the value of e = 0.866. It can also be shown that in that 
case the deformation parameter would actually be e 2 , so 
that for the same deformation we would have e 2 = 0.75 
which is obviously much larger than the parameter which 



we have actually used here. Such a high value of the 
deformation parameter goes against the very essence of 
the perturbative nature of the method. This means that 
while we have terminated the Fourier series and also the 
eigenvalues at the second order of smallness when work- 
ing with A = ^pj, for A = e 2 we would have to consider 
higher order terms to get the same accuracy. Finally, 
we note that the same formalism can also be adopted 
by perturbing the boundary of a square or a rectangular 
membrane for which the results are exactly known. 
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FIG. 1: Comparison of the energy eigenvalues obtained nu- 
merically and analytically for a particle enclosed in a super- 
circular boundary (in units of 2t ^ fl a ) 
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